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Abstract 

The transport cross-section based on inflow transport approximation can significantly improve the 
accuracy of light water reactor (LWR) analysis, especially for the treatment of the anisotropic 
scattering effect. The previous inflow transport approximation is based on the moderator cross- 
section and normalized fission source, which is approximated using transport theory. Although the 
accuracy of reactivity is increased, the Po flux moment has a large error in the Monte Carlo code. 
In this study, an improved inflow transport approximation was introduced with homogenization 
techniques, applying the homogenized cross-section and accurate fission source. The numerical 
results indicated that the improved inflow transport approximation can increase the Po flux 
moment accuracy and maintain the reactivity calculation precision with the previous inflow 
transport approximation in typical LWR cases. In addition to this investigation, the improved 
inflow transport approximation is related to the temperature factors. The improved inflow 
transport approximation is flexible and accurate in the treatment of the anisotropic scattering 
effect, which can be directly used in the temperature-dependent nuclear data library. 

Keywords: Inflow transport approximation; Anisotropic scattering effect; Homogenization 
techniques; Light water reactor 


1. Introduction 

In light water reactor (LWR) analysis, the anisotropic scattering effect is generally treated by a 
simplified approximation, such as the transport cross-section (XS) with isotropic scattering 
approximation from the computational efficiency viewpoint. This approximation has been used in 
general LWR problems and has reached operational experiences, regardless of its simplicity. 
However, recent studies suggest that the conventional approximation is inaccurate due to several 
problems, such as high-leakage, small cores, and fuel assemblies with control rods [1][2]. 
Therefore, the treatment of the anisotropic scattering effect is important for the recent lattice 
physics code. 

There are three major issues in the treatment of the anisotropic scattering effect. The first is a 
higher-order scattering matrix. The P» or higher-order scattering matrix is sufficient to maintain 
calculation accuracy [3]. However, a large amount of memory and a long computing time are 
required to incorporate the higher-order scattering matrix. The high-order scattering matrix is 
implemented for the whole-core calculation in the Michigan Parallel Characteristics Transport 
(MPACT) code [4], and whole-core problems are calculated using the P, scattering approximation. 
The second is the application of the angular-dependent total cross-section. Multi-group total XS is 
obtained from angular flux condensation [5][6][7]. Conversely, it employs the time-consuming 
Monte Carlo code, and it is also very difficult to attain the accurate angular flux. A previous study 
showed that 160 batches with 1 billion neutrons per tally can reach the converged results [8], 
which is also an impractical method. The third is the transport approximation for transport XS [9]. 
It can be used in the Po scattering matrix and is efficient for large-core calculations. Nevertheless, 
the key point in the inflow transport approximation is the accurate higher-order flux moment from 
the neutron transport calculation. Yamamoto [2] adopted the typical neutron spectrum for the Po 
flux moment, and the higher order flux moment was inversely proportional to the total XS but did 
not improve the accuracy of reactivity compared to that using the outflow transport 
approximation. Choi [10] developed an inflow transport approximation from a one-dimensional 
Py transport equation, improving the accuracy in cases of high anisotropic scattering. However, 
Choi’s one-dimensional Py transport equation is based on the moderator region and the 
normalization of the fission source. Although the accuracy of the reactivity calculation in the high 
anisotropic scattering cases improved, it was approximated with the theoretical model. 


Nevertheless, there are many other studies related to the neutron transport problem and cross- 
section [11][12][13]. 

In this study, the improved inflow transport approximation was investigated, and its 
application in different LWR cases was tested. The transport XS was computed by solving the 
one-dimensional Py transport equation with the homogenized cross-section and accurate fission 
source, verified through typical LWR cases. The accuracy and application of these transport 
approximation methods can be determined by comparing the k-infinity, k-effective, and Po flux 
moments to the reference results. In addition to the inflow transport approximation investigation, 
the application of the improved inflow transport approximation was tested using different typical 
LWR cases. Moreover, the results of inflow transport approximations were performed in the lattice 
code DRAGONS.0.5 [14] and compared to the results of experimental values and Monte Carlo 
code cosRMC [15]. The numerical results showed that the improved inflow transport 
approximation can reach a higher accuracy than the previous method for the Py flux moment, and 
can retain the same accuracy as the previous inflow transport approximation for the reactivity 
calculation. In addition, the improved inflow transport approximation is only related to the 
temperature factor through the verification of typical LWR cases. As a result, the improved inflow 
transport approximation is precise and flexible in the LWR analysis when applied to a 
temperature-dependent nuclear data library. 

The remainder of this manuscript is organized as follows. Section 2 describes the theory and 
methodology of the different inflow transport approximations. Section 3 explains the discrepancy 
analysis between these two methods. Section 4 presents the numerical results for the different 
LWR cases. The conclusions are expressed in Sect. 5. 


2. Theory and methodology 
2.1 Previous inflow transport approximation 

The previous inflow transport approximation was developed by Choi, which is based on the 
moderator region and normalized fission source. In this method, the multi-group form of the one- 
dimensional Py transport equation is derived as follows: 


(0) 
Where 
(0) 
(0) 
(0) 
(0) 
where is the angular flux in group g, is the cosine angle in the z-direction, is the Legendre 
polynomial, / is the order of the Legendre polynomial, is the /-th moment macroscopic XS of x 
reaction (f: fission, t: total, s: scattering), and is the /-th flux moment. 
It should be noted that the distinction between and (n=1,2...) disappears in the absence of 
self-shielding (i.e., infinity dilution) [16]. 
The multigroup transport equation in Eq. (1) can be extended using Eqs. (2)-(5), as follows: 
. (0) 
Multiplying Eq. (6) by the Legendre polynomial and integratingfrom - 1 to +1, and using the 
following Bonnet’s recursion formula and the orthogonality of the Legendre polynomials: 


; (0) 

(0) 
Inserting Eqs. (7) and (8) into Eq. (6) will lead to the following equation: 

(0) 
Approximating the spatial shape using buckling B: 

; (0) 
Inserting Eq. (10) into Eq. (9) will produce the following equation: 

(0) 


Additionally, another approximation was introduced for Choi’s method. The fission source is 
forced to normalize, and the value of k is set to 1, leading to: 


(0) 
(0) 


In Choi’s method, Eq. (11) is then converted to the following equation: 


Solving Eq. (13), the then leads to the following: 


(0) 
Where the transport XS is written as: 
- (0) 
When n is an even number and n-1 and n+1 are odd numbers. 
For n-1: 
(0) 
For n+1: 
l (0) 
Inserting Eqs. (16) and (17) into Eq. (14) leads to the following: 
(0) 
Where 
(0) 
(0) 
(0) 


If the order of n is the maximum value, becomes 0. For the odd n, Eq. (14) becomes the 
following: 
; (0) 
The high-order flux moments and Po transport XS are then generated. Furthermore, the Po 
scattering matrix is calculated using the transport XS, and can be used in the transport 
computation. Only the corrected diagonal elements are as follows: 
: (0) 
The transport equation is related to the geometric value z. It was calculated for a dominated 
moderator material in each problem case, and showed a problem-dependent result. 
2.2 Improved inflow transport approximation 
From the above derivation, a key point in the inflow transport approximation is applying the 
accurate flux moment from a one-dimensional Py transport equation, considering the leakage 
effect and anisotropic scattering effect. This is the same as the assembly homogenization 
techniques [17]. However, there are two terms (multi-group data and fission source) that are 
inconsistent with the theoretical model through the derivation of the previous inflow transport 
approximation. 
When applying the flux shape with buckling B, the multigroup XS and scattering matrix in 
Eqs. (3) and (4) become the following: 
(0) 
(0) 
This shows that applying the flux shape with buckling B eliminates the geometrical value z in 
the derivation. Therefore, the previous one-dimensional Py transport equation becomes a 
homogenized system. The cross-section and scattering matrix data are homogenized in the entire 
region, such as fuel, cladding, and moderator regions. 
Excluding the homogenized system, the accurate fission source is also treated using 
homogenization techniques. Therefore, the transport equation in Eq. (11) is replaced by the 
following: 


(0) 
The solution of flux moment then leads to the following: 
(0) 
Where 
(0) 
(0) 
When n is an even number and n-1 and n+1 are odd numbers. 
For n-1: 
(0) 
For n+1: 
(0) 
Inserting Eqs. (30) and (31) into Eq. (27) leads to the following: 
(0) 
Where 
(0) 


(0) 
(0) 
If the order of n is the maximum value, becomes 0. For the odd n, Eq. (27) becomes the 
following: 
; (0) 
Therefore, the flux moment and transport XS were obtained using the improved inflow 
transport approximation. Compared with the form of the previous inflow transport equation, the 
improved inflow transport equation was based on the homogenization XS and accurate fission 
source, which is consistent with the homogenization techniques. 
2.3 Implementation of improved inflow transport approximation 
In Choi’s parametric study, the inflow transport approximation is insensitive to the buckling 
and order of the Legendre polynomial. For the value of buckling, even with a ten-fold difference 
in buckling, the influence of k-effective is less than 2 pcm error. Therefore, the typical buckling 
value (B’=0.0001) was used for the calculation. For the order of the Legendre polynomial, the P: 
transport equation is sufficient to calculate the converged Po transport XS, which has the same 
impact as the P; transport equation [9]. 
In the improved inflow transport approximation, the typical buckling value and P; transport 
equation were used. The detailed iteration process is as follows. 
(a) Obtain the homogenization data based on the outflow transport approximation in an actual 
problem. The DRAGON code produces the following homogenization data: 
re (0) 
s> (n=1,2) (0) 
where 7 is used for every region, mod is the moderator region, x is represented as the reaction (f: 
fission, t: total, s: scattering), and n is the order of the Legendre polynomial (the anisotropic 
scattering effect for non-light nuclides is ignored). 
(b) Initial flux moments and transport XS. 


, (n=1,2), (n=0,1,2) (0) 
(c) Update the Py and P, flux moments using Eq. (32), and the flux moments converge. 
(n=0 and 2) (0) 
(d) Update the P, flux moments. 
(0) 
(e) Update the Py transport XS. 
(0) 


(£) Repeat steps (b)-(e) until the flux moments and transport XS converge. 

In the above procedure, the flux moments and transport XS may be negative. The iterative 
equation is as follows:. 

(0) 
(0) 
where k is the iterative number. 

In this study, the investigation focused only on the dominant moderator material (H20). 
Therefore, the transport XS for H:O was selected for the improved transport approximation. A 
detailed discrepancy analysis between these transport approximations is presented in the next 
section. 


3. Discrepancy analysis 

The typical LWR benchmark was selected in this study to determine the accuracy of these 
transport approximations. This is the same as the problem Bettis Atomic Power Laboratory 
(BAPL-1) in the WIMS-D Library Update Project (WLUP) [18]. The geometry is shown in Fig. 1, 
and the material composition is listed in Table 1. The temperature was set to 300 K for all regions, 
and the inter-gap was omitted. 


Fig. 1 The geometry of BAPL-1 


Table 1 The material composition of BAPL-1 


Material Nuclide Density (atom/barn.cm) 
Moderator 'H 6.6760x10° 

1O 3.3380x10? 
Cladding Al 4.8994x10° 


Fuel SU 3.1120x104 
238 2.312710" 
16O 4.6946x10° 


Figure 2a shows the transport XS for H20 in three transport approximations: outflow, previous 
inflow, and improved inflow transport approximation. The outflow transport approximation had a 
large discrepancy in the thermal and fast energy ranges to the inflow transport approximation. The 
transport XS for H:O showed a similar tendency between these two inflow transport 
approximation methods. 


(a) (b) 
(c) 


Fig. 2 (Color online) Comparison of transport XS and normalized Po flux moment in the BAPL-1 case. (a) 
Transport XS for H20 ; (b) Comparison of normalized Po flux moments in moderator region; (c) Comparison of 
normalized Po flux moments in the homogenization region 


Furthermore, to determine the difference between these inflow transport approximations, a 
comparison of the Po flux moment is shown in Figs. 2b and 2c. The treatment of the fission source 
and XS were the main differences between these inflow transport approximations. The reference 
result was obtained using the Monte Carlo code cosRMC, which was developed by Tsinghua 
University and the State Power Investment Corporation Research Institute (SPIC). Therefore, 
there are two types of Po flux moment comparisons in this study. 

1. The comparison in the moderator region. 

(a) Normalized fission source + moderator XS (previous method). 

(b) Accurate fission source + moderator XS. 

(Reference results: cosRMC’s Po flux moment in the moderator region) 
2. The comparison in the homogenized region. 

(a) Normalized fission source + homogenized XS. 

(b) Accurate fission source + homogenized XS (improved method). 
(Reference results: cosRMC Po flux moment in the homogenized region) 

In Fig. 2b, the normalized Ps flux moment in the previous inflow transport approximation is 
significantly different from the cosRMC code. For the moderator system, the P, flux moment 
changed with different fission sources. When adopting the normalized fission source, the Po flux 
moment became more difficult than the accurate fission source. In Fig. 2c, the normalized Po flux 
moment in the improved inflow transport approximation agrees well with the cosRMC code. For 
the homogenized system, the Po flux moment also changed with different fission sources, which 
had the same tendency as the moderator system. Moreover, for the same fission source with 
different XS cases, the homogenized XS lead to a 1 MeV decrease in the Py flux moment and was 
closer to the LWR system [19], but the P, flux moment in the moderator XS was different. 
Therefore, the inflow transport approximation, with an accurate fission source and homogenized 
XS, can improve the accuracy of the Po flux moment compared with that using the previous 
method. 


4. Numerical results 
4.1 General description of verification 

The lattice code DRAGONS.0.5 was used for verification. It was developed by Ecole 
Polytechnique de Montréal for lattice calculations, solving the neutron transport equation with a 
deterministic approach, and can read different types of multi-group libraries, such as the versions 
of WIMS-D and MATXS. In this study, the WIMS-D 70 library [20] based on ENDF/B-VII.0 [21] 
was selected to study the influence of the anisotropic scattering effect. Different transport 
approximations are implemented in the WIMS-D 70 library through the transport XS and 
scattering matrix. Resonance self-shielding is an equivalence theory (SHI+LJ+LEVEL2) with the 
generalized Stamm’ler method [22], Livolant and Jeanpierre normalization factor, and Riemann 
integration [23]. The neutron transport calculation method was the interface current method 
(SYBILT). cosRMC was used as the reference code. The results of the A-infinity, k-effective, and 
P, flux moments are compared for various verification problems. 
4.2 Typical LWR fuel cell and fuel assembly cases 


A typical LWR fuel cell and fuel assembly were used to verify the application of different 
transport approximations. The factors of the fuel cell, fuel assembly, boron concentration, fuel 
enrichment, control rod material, and temperature in the fuel assembly with a control rod were 
selected as verification cases. 


Fig. 3 (Color online) Geometry shape of 1 A(left) and 2A (right) 

The single fuel cell and 17x17 fuel assembly benchmark were extracted from problem-1A and 
2A in VERA [24], which contained 3.1 wt% UO; fuel, 0.743 g/cc moderators with a 1,300 ppm 
boron concentration, and Zr-4 for cladding. The temperature was 300 K for all regions. The 
geometric shapes for 1A and 2A are shown in Fig. 3, and a detailed description of each case is 
given in Table 2. 

Table 3 shows the kint results for each case with different transport approximations. Compared 
with the kint in cosRMC, the kint error in DRAGON for inflow transport approximations was less 
than 100 pem. Contrarily, the error of king in DRAGON for outflow transport approximation was 
less than 100 pcm in the fuel cell and fuel assembly, but could cause approximately 150 to 250 
pcem in fuel assembly with the control rod. This discrepancy is consistent with that of previous 
research [9]. 

Table 2 Detailed description of each case 


Cases Description 2U w/o Temperature (K) Moderator 

Fuel Clad Mod g/cm? ppm 
1A fuel cell 3.1 300 300 300 0.743 1300 
2A fuel assembly 3.1 300 300 300 0.743 1300 
2H-typical fuel assembly+24 B,C 3.1 300 300 300 0.743 1300 

2H-0 ppm 2H-tyical+Oppm 3.1 300 300 300 0.743 0 
2H-4.8% 2H-tyical+4.8 % 4.8 300 300 300 0.743 1300 
2G fuel assembly+24 AIC 3.1 300 300 300 0.743 1300 
2H-600 K 2H-tyical+600 K 3.1 600 600 600 0.743 1300 


Compared with the knr in the previous inflow transport approximation, the results of the 
improved transport approximation were less than 10 pcm in cases 1A and 2A. Moreover, the 
results of the improved transport approximation were underestimated by 11 to 49 pcm for 2H- 
typical, 2H-0 ppm, 2H-4.8%, and 2G, and was also overestimated by 46 pem for 2H-600 K. This 
tendency indicates that the transport approximations have less influence on the fuel cell and fuel 
assembly, but they have some effect on the fuel assembly with control rod cases. For the 300 K 
cases, the application of homogenization techniques in the improved transport approximation is 
based on the homogenous XS and accurate fission source, considering the influence of the 
resonance self-shielding effect in fuel materials. Compared with the previous inflow transport 
approximation, the improved transport approximation can absorb neutrons in different cases, 
which finally decreases the kint results. During the 600 K case, the impact of the resonance self- 
shielding effect was weakened in the improved inflow transport approximation with an increase in 
temperature [25], leading to an increase in kins compared to that in the previous transport 
approximation. 

Table 3 Results of Aine in different cases 


Case cosRMC DRAGON 

Method! Method2 Method3 

1A 1.20116 1.20037 1.20033 1.20038 
(+30 pem) (-80 pem) (-83 pem) (-78 pem) 

2A 1.19391 1.19404 1.19405 1.19407 
(+30 pem) (13 pem) (14 pem) (16 pem) 

2H-typical 0.81444 0.81571 0.81422 0.81393 
(+31 pem) (127 pem) (-22 pem) (-51 pem) 

2H-0 ppm 0.89927 0.90174 0.89999 0.89973 
(+32pcem) (247 pem) (72 pem) (46 pem) 

2H-4.8 % 1.00206 1.00419 1.00260 1.00211 

(+31 pem) (213 pem) (54 pem) (5 pem) 

2G 0.87892 0.88044 0.87893 0.87882 
(+29 pcm) (152 pem) (1 pem) (-10 pem) 


2H-600 K 0.78723 0.78956 0.78720 0.78766 

(427 pem (233 pem) (-3 pem) 43 pcm 
*Method1: outflow transport approximation; Method2: previous inflow transport approximation; Method3: 
improved inflow transport approximation 


(a) (b) 

(c) (d) 

(e) (f) 
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Fig. 4 (Color online) Comparison of normalized Po flux moment with the previous method. (a) 1A fuel cell; (b) 
2A fuel assembly; (c) 2H-typical fuel assembly with B,C; (d) 2H-0 ppm fuel assembly with B,C at 0 ppm; (e) 2H- 
4.8% fuel assembly with B,C in 4.8% fuel enrichment; (f) 2G fuel assembly with AIC; (g) 2H-600 K fuel 
assembly with B,C in 600 K 


As a result, the different inflow transport approximations reached a higher accuracy than the 
outflow transport approximation in reactivity calculations, and the error of kint was less than 50 
pem for different inflow transport approximations. 

Figures 4 and 5 show the comparison of the normalized Py flux moment with different inflow 
transport approximations for these cases. Specifically, compared with the results of cosRMC, the 
Po flux moment in the improved inflow transport approximation had a better agreement than the 
previous method. Furthermore, the Po flux moment of the previous inflow transport approximation 
was slightly overestimated by approximately 1 MeV and underestimated in the thermal energy 
range. Particularly, the previous inflow transport approximation made the Po flux moment more 
difficult than the improved method. 

In addition to the above investigation, the results of the improved inflow transport 
approximation in different LWR cases were in good agreement with the cosRMC, which is more 
versatile than the previous method. It can achieve higher accuracy for the cases of the fuel cell, 
fuel assembly, different boron concentrations, fuel enrichment, control rod material, and 
temperature in fuel assembly with the control rod. Therefore, the improved inflow transport 
approximation is accurate, universal, and suitable for LWR analysis. 


(a) (b) 
(c) (d) 
(e) (f) 

(g) 


Fig. 5 (Color online) Comparison of normalized Po flux moment in the improved method. (a) 1A fuel cell; (b) 
2A fuel assembly; (c) 2H-typical fuel assembly with B,C; (d) 2H-0 ppm fuel assembly with B,C at 0 ppm; (e) 2H- 
4.8% fuel assembly with B4C in 4.8% fuel enrichment; (f) 2G fuel assembly with AIC; (g) 2H-600 K fuel assembly 
with B.C in 600 K. 


4.3 Analysis of influence factors 

In a previous investigation, the inflow transport approximation was problem-dependent. 
However, transport approximations are problem-independent data in most nuclear data libraries. 
Hence, we need to study the influence of the transport approximation and produce a nuclear data 
library with a general transport approximation. 


(a) 
(b) 


Fig. 6 (Color online) Macroscopic transport XS (left) and relative error (right) of 2A for H:O. (a) Results of 
previous inflow transport approximation; (b) Results of the improved inflow transport approximation. 


Figure 6 shows the macroscopic transport XS and relative error of 2A for H20 in the different 
inflow transport approximations. The macroscopic transport XS for H2O was similar in most LWR 
cases, but the biggest difference existed in the thermal energy group in the 600 K case. The 
relative error of transport XS in the 600 K case was approximately 30% to 40% in the 2A case in 
the different inflow transport methods, but was less than 20% in other cases. Therefore, the most 
significant factor influencing inflow transport approximations is temperature through the 
comparison of transport XS. 

Moreover, for the previous inflow transport approximation, the transport XS was based on the 
moderator region, and its relative error to 2A was less than 1% with the same moderator XS. 
Therefore, the previous inflow transport approximation was only related to the moderator region. 
For the improved inflow transport approximation, the transport XS was based on the 
homogeneous system, and its relative error to 2A was over 5% in 0.28 eV to 27.7 eV. Specifically, 
in a homogeneous system, the resonance self-shielding effect for fuel and absorber nuclides was 
considered in the derivation of transport XS. Then, the neutron spectrum was changed according 
to the different nuclide density ratios in the fuel and absorber nuclides. Thus, the transport XS in 
the improved inflow transport approximation was different from the previous inflow method, 
which is consistent with the actual LWR model. 

The transport XS had a strong correlation with temperature and had less difference in the 
different LWR operation factors, such as fuel assembly, boron concentration, fuel enrichment, and 
control rod materials in the fuel assembly with the control rod. There are two different WIMS-D 
70-group libraries, based on the improved inflow transport approximation, to test the influencing 
factors of transport approximations. The first is a problem-dependent library. The transport XS and 
scattering matrix data for 'H and "O are problem-dependent, which are developed by the 
improved inflow method in actual cases. The second is a problem-independent library. The 
transport XS and scattering matrix data for H-1 and O-16 are problem independent, developed by 
the improved inflow method, from the typical 2H case (BuC, 300 K, 1,300 ppm, and 3.1 wt%). 

Table 4 shows the results of the kint in the improved transport approximation. Compared to 
the kine results between the problem-dependent and problem-independent libraries, the kint error 
was less than 20 pem for most cases. However, the kint error was 37 pcm for the 2H case at 600 K. 
This indicates that the inflow transport approximation considering fuel and absorber nuclides has 
less influence on kint calculations, but the factor of temperature can influence the kint results. This is 
consistent with the difference in the transport XS in Fig. 6. Although the largest difference of 
transport XS was in the resonance energy range, the magnitude of transport XS in the resonance 
energy range was comparatively smaller than the thermal energy range, and its influence on the kint 
error was less than 20 pcm in each case. Therefore, the improved inflow transport approximation 
is only related to the temperature factor, and it can be easily implemented in a temperature- 
dependent nuclear data library. 

Table 4 Results of kine in improved inflow transport approximation 


Cases cosRMC DRAGONS.0.5 
Dependent Independent 

2H(typical) 0.81444 0.81393 0.81393 
(+31 pcem) (-51 pem) (-51 pem) 

1A 1.20116 1.20038 1.20040 
(+30 pem) (-78 pem) (-76 pem) 

2A 1.19391 1.19407 1.19408 
(+30 pem) (16 pem) (17 pem) 

2H(0 ppm) 0.89927 0.89973 0.89967 
(+32 pem) (46 pcm) (40 pcm) 

2H(4.8 wt%) 1.00206 1.00211 1.00231 
(+31 pem) (5 pem) (25 pem) 

2G(AIC) 0.87892 0.87882 0.87880 
(+29 pem) (-10 pem) (-12 pem) 

2H(600K) 0.78723 078766 0.78729 

(+27 pem) (43 pem) (6 pem) 


4.4 Application of inflow transport approximation to core cases 

The Babcock & Wilcox (B&W) critical experiment benchmarks [26] selected to verify the 
application of different transport approximations in core cases. The B&W 1484 core 1 and core 3 
cases are selected in this section, and the geometry shape is given in Fig. 7. The B&W 1484 core 1 
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is the 48x48 case, and the B&W 1484 core 3 is the 68x68 case. 

Table 5 shows the results of the B&W 1484 core 1 and core 3 cases (2D and 3D). The results 
indicated that the outflow transport approximation can introduce an error of over 1,000 pem in 
these cases, and both inflow transport methods were close to the reference values. For the B&W 
1484 core 1 (2D and 3D) cases, the results of kint in the improved inflow transport approximation 
were overestimated by approximately 364 pcm and 291 pcm, respectively. For the B&W 1484 
core 3 (2D and 3D) case, the results of kint in the improved inflow transport approximation were 
overestimated by approximately 115 pcm and 121 pcm, respectively. The knr error in the B&W 
1484 core 1 and core 3 cases was inversely proportional to the cases of the LWR assembly. In 
addition, the different transport approximations were sensitive to the small core cases, which can 
cause over 1,000 pcm for the knr calculations in core 1 and approximately 600 pcm in core 3. 

Although the previous inflow method is close to the reference results, the improved inflow 
method is close to the actual system, which is consistent with the theoretical model. The reason for 
this tendency may come from the limitation of the WIMS-D library, such as ignorance of 
scattering resonance integral, lumped fission spectrum, etc. 

In summary, the outflow transport approximation can introduce an error of over 1,000 pem in 
the B&W 1484 core 1 and core 3 cases, but the inflow transport methods are relatively accurate in 
the core cases. 


Fig. 7 (Color online) Geometry shape of B&W 1484 Core 1 (left) and Core 3 (right) 


Table 5 Results of Kinser in B&W 1484 Core 1 and Core 3 


Case CosRMC/ DRAGON 

Experiment Method1 Method2 Method3 

Core 1(2D) 1.01270 1.02851 1.01198 1.01562 
(+30 pem) (1581 pem) (-72 pem) (292 pem) 

Core 1(3D) 1.00020 1.01387 0.99765 1.00076 
(-) (1367 pem) (-235 pem) (56 pem) 

Core 3(2D) 1.01630 1.02800 1.02183 1.02298 
(+30 pem) (1170 pem) (553 pem) (668 pem) 

Core 3(3D) 1.00000 1.01331 1.00725 1.00846 
- 1331 pem 725 pem 846 pem 


*Method1: outflow transport approximation; Method2: previous inflow transport approximation; Method3: 
improved inflow transport approximation 


5. Conclusion 

In this study, the improvement and application of inflow transport approximation in LWR 
analysis were performed by comparing the results of the Po flux moment, ker, and kine with the 
reference results. The improved inflow transport approximation was based on the homogenized 
cross-section and accurate fission source, which was different from the previous method. 

For the LWR assembly cases, the results of three transport approximations (outflow, previous 
inflow, and improved inflow) were compared for the cases of the fuel cell, fuel assembly, boron 
concentration, fuel enrichment, control rod materials, and temperature in the fuel assembly with a 
control rod. The Kine error in the inflow transport approximation was less than 100 pem compared 
to that in the Monte Carlo code cosRMC, which was more accurate than the outflow transport 
approximation in the treatment of the anisotropic scattering effect. Compared to the Po flux 


moment, the results of the improved inflow transport approximation were closer to the cosRMC 
code in each case. The homogenized cross-section and accurate fission source were the main 
factors influencing the Po flux moment. In addition to the above investigation, the improved 
inflow transport approximation is only related to the temperature factor for typical LWR cases. 
Excluding the temperature case, the Aine error was less than 20 pcm between the problem- 
dependent and problem-independent transport approximations. For the B&W 1484 core cases, the 
outflow transport approximation introduced an error of over 1,000 pcm. The inflow transport 
approximation significantly improved the kint accuracy, and the results were closer to the reference 
value. 

In summary, the improved inflow transport approximation is more flexible and accurate in the 
treatment of the anisotropic scattering effect, and it can be easily implemented in the temperature- 
dependent nuclear data library for LWR analysis. 
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